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Most experiments that search for direct interactions of WIMP dark matter with a target can 
distinguish the dominant electron recoil background from the nuclear recoil signal, based on some 
discrimination parameter. An acceptance region is defined in the parameter space spanned by the 
recoil energy and this discrimination parameter. In the absence of a clear signal in this region, a 
limit is calculated on the dark matter scattering cross section. Here, an algorithm is presented that 
allows to define the acceptance region a priori such that the experiment has the best sensitivity. 
This is achieved through optimized acceptance regions for each WIMP model and WIMP mass that 
is to be probed. Using recent data from the CRESST-II experiment as an example, it is shown that 
resulting limits can be substantially stronger than those from a conventional acceptance region. In 
an experiment with a segmented target, the algorithm developed here can yield different acceptance 
regions for the individual subdetectors. Hence, it is shown how to combine the data consistently 
within the usual Maximum Gap or Optimum Interval framework. 



I. INTRODUCTION 

Although the existence of dark matter is now well es- 
tablished (see e.g. [1]), we have not yet succeeded in de- 
termining its nature. A highly motivated class of models 
predicts dark matter to be in the form of Weakly Inter- 
acting Massive Particles (WIMPs) (see e.g. [2]). Direct 
searches for scatterings of such WIMP dark matter par- 
ticles off nuclei (see e.g. [3]) probe and constrain models 
in highly relevant regions of parameter space. Limits on 
the cross section of the scattering process are calculated 
based on the exposure of the experiment as well as events 
that are observed in an acceptance region. 

The question arises how this acceptance region should 
be defined. Here, an algorithm is presented that allows 
to define an optimum acceptance region such that the ex- 
periment has the best sensitivity for a given WIMP sig- 
nal. Hence, one can expect to obtain the most stringent 
limit on the dark matter scattering cross section. The 
algorithm is illustrated using recent data from CRESST- 
II [4], but can be employed in any experiment. It is 
also shown how to consistently combine data from in- 
dividual detectors or experiments with differing spectral 
acceptance in the framework of Yellin's Maximum Gap 
or Optimum Interval methods [5] . 



II. THE OPTIMUM ACCEPTANCE REGION 

Figure 1 shows data collected in the CRESST-II ex- 
periment [4] that will be used as an example to explain 
the method. For each particle interaction in the target, 
the experiment collects two signals: The recoil energy E 
is taken from the phonon detector, and the light output 
L is measured by a separate light detector in units of 
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keV electron equivalent keV co , defined such that 122 keV 
gammas from a 57 Co calibration deposit 122 keV ee in the 
light detector [6] . The light yield y is defined as the ratio 
of these two parameters, y = L/E. It serves as discrim- 
ination parameter to distinguish the dominant electron 
recoil background in the band around y « 1 keV cc /keV 
(figure 1) from the nuclear recoil signal which is expected 
around zero light yield. Other experiments have differ- 
ent discrimination parameters such as the charge yield 
Q/E [7, 8] or the ratio of delayed and prompt scintilla- 
tion S^/Si [9, 10]. The algorithm presented in the fol- 
lowing is applicable irrespectively of the nature of the 
discrimination parameter. 

Of course any algorithm that aims to define an ac- 
ceptance region needs to be blind to the particular dis- 
tribution of individual events in the data. To facilitate 
this, the data can be regarded as a background density 
Qh(y,E) which is calculated based on a parametrization 
of the data. Hence, g^(y,E) is a real valued function 
defined for the whole parameter space (y, E) and is inde- 
pendent of individual signal candidate events. 

In the CRESST-II example considered here, the back- 
ground can be modeled as a Gaussian band, prominent 
in figure 1. The energy dependence of the background 
is taken from the measured spectrum [11]. To find the 
mean L(E) and width a(E) of the background band, 
the data is fitted with a Gaussian function that allows 
for the observed scintillator non-proportionality in the 
mean [12], and with a parametrization for the width 
a 2 (E) = o-Q+alE+<7%E 2 as used by the collaboration [6]. 

In addition to Qb, there is an expected signal density 
(JQ s (y, E) that is proportional to the WIMP-nucleon scat- 
tering cross section a. This density is calculated from the 
position and width of the nuclear recoil band [4] and the 
expected WIMP spectrum oWj&E [13, 14]. 

Only events that are found within the acceptance re- 
gion are considered to calculate a limit. In the CRESST- 
II example, the acceptance region is parametrized by its 
upper boundary y mli x{E), since the negative boundary 
can be taken to be at y m in(E) = const = — lkeV 00 /keV 
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FIG. 1: Two-dimensional histogram of events collected by one 
CRESST-II detector module in the light yield-energy-plane 
y(E) (detector Verena after an exposure of 24.11 kg d [4]). 
Counts per bin are color coded according to the grey scale 
on the right. The band around a light yield of unity comes 
from the dominating electron and gamma background. The 
dark matter induced nuclear recoil signal is expected around 
zero light yield. Entries with negative light yield originate 
in the amplitude fitting procedure that allows for negative 
amplitudes in order to treat noise in an unbiased way [6] . The 
hatched area (green) is an example of the optimum acceptance 
region for a 100 GeV/c 2 WIMP, and the area at lower energies 
(blue) an example for the case of a 10 GeV/c 2 WIMP. 



which is equivalent to y m i n (i?) = — oo for all practical 
purposes. For a particular realization of the acceptance 
region {y meLX (E), y m - ln (E)}, the expectation values for the 
observed number of signal and background events are 
then simple integrals 

dE / dy <TQ s (y,E) (1) 

■://,,; = I dE dy g h (y,E). (2) 

Jy m in{E) 

Given the acceptance region {y m ax{E),y m in(E)}, a 
particular experiment will observe (after unblinding) a 
certain number of accepted events n b s £ N. This num- 
ber then allows to calculate some upper limit n Sj go still 
compatible with the data at a stated confidence, typically 
90%, often using the Maximum Gap or Optimum Inter- 
val methods (section III) . Here it suffices to note that the 
limit depends on the acceptance region: 

« s ,9o = «-s,9o(«obs) = n St go(y max (E),y m ia{E)). (3) 

To transfer this number into a limit ago on the cross sec- 
tion, equation 1 is evaluated for an arbitrary a', yielding 
(rig), and 090 can then be calculated according to 



^90 
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n s,90 



Since (n' s ) oc a', the primed variables eventually drop out 
of the equation again. 

The upper limit on the cross section ago will of course 
depend on the particular realization of the events in a 
given experiment, so it cannot be used as the objective 
function. However, the expectation value (ergo) (the sen- 
sitivity of the experiment) is independent of a given ex- 
perimental outcome and can be used instead. Its value 
follows from equation 4: 
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The upper limit n Sj go depends on the number of ob- 
served events n b s , which in turn depends on the number 
of background events rib- The distribution of rib can be 
assumed to be Poissonian with expectation value (rib) 
from equation 2. To calculate the sensitivity (ergo), we 
can then simply use the definition of the expectation 
value: 

, 00 

(ego) = X! n s,9o(™obs) P(n ohs ) 
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The sensitivity (ergo) depends on the acceptance region 
{yma.x(E), y m i n (E)} through equations 1 and 2 and is the 
objective function of choice. 

An algorithm for finding the optimum acceptance re- 
gion varies y ma , x (E) and y m i n (E) until (ago) is minimal. 
To this end, the (y, E) plane is binned and varied on its 
borders: For each energy bin, the acceptance region is in- 
creased or decreased as long as it improves the sensitivity. 
The algorithm is iterated until the sensitivity converges. 

Some technical remarks will help to implement the al- 
gorithm. Although e _ ( nb ' can be drawn out of the sum 
in equation 6, this does not mean that (ago) — > as 
(rib) - * 00, since n s ,90 counteracts. Instead, the sum 
is very well behaved, as can be seen in figure 2, where 
n s go(n bs) is calculated using Poisson statistics. Had 
one precise knowledge of the background, n Sj go could also 
be calculated e.g. from the Fcldman-Cousins scheme [15] 
(dashed line in figure 2), but this is not used here. Since 
the sum is well behaved, it can be tabulated for integer 
(rib) and then interpolated to speed up its computation. 
Given the simple Gaussian behavior of the signal and 
background densities g s and gy,, it is not surprising that 
the objective function is also well behaved, so the algo- 
rithm rapidly converges to the optimum acceptance re- 
gion. Eventually, this region can be decreased again by a 
small amount (e.g. by 0.1% of the accepted WIMP spec- 
trum) to prevent numerical ambiguities from showing up 
at higher energies. The optimum acceptance region then 
naturally extends to a maximum energy above which the 
acceptance is zero, y max = y min . 

Since the acceptance region is calculated using the sig- 
nal expectation ag s (y,E), it will be different for each 
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FIG. 2: The expectation value (n a ,go). The solid curve gives 
the employed curve based on Poisson statistics as shown in 
equation 6; dashed the (not used) variant in which the back- 
ground is taken into account using the ordering scheme of 
Feldman and Cousins. 



WIMP model that is probed. In particular, a separate 
acceptance region is calculated for each WIMP mass. 
For the CRESST-II example considered here, two opti- 
mum acceptance regions are shown in figure 1 for generic 
WIMPs with masses of 10GeV/c 2 and 100GeV/c 2 . For 
the latter, the acceptance region follows the resolution 
of the light channel. The acceptance regions are not 
smoothly bounded but show small dents. This is because 
the background density Qb needed to evaluate equation 2 
is calculated using the observed spectrum, smeared in y 
by the resolution of the light channel. Hence, the accep- 
tance region is reduced for example around 11 keV where 
a gamma line appears in the spectrum [11]. 

It may be surprising to note the shape of the opti- 
mum acceptance region for low mass WIMPs: It is in 
fact beneficial to increase the acceptance toward lower 
recoil energies. This is due to the recoil spectrum being 
a sharply falling exponential, confined to the lowest ener- 
gies. Hence, the increasing electron/gamma background 
in this parameter region is being overcompensated for. 

The improvement brought by this algorithm can be il- 
lustrated calculating a limit based on the above data set. 
To this end, an analysis threshold of 3keV is imposed 
to stay well above detection threshold [16], and the same 
WIMP expectation as used by the CRESST collaboration 
is probed [6]. In particular, the WIMPs are expected to 
be distributed in the Milky Way in an isothermal halo 
with a local density of 0.3 GeV/c 2 /cm 3 , through which 
Earth moves with a velocity of 230km/s [13]. Recoils 
with high momentum transfers are suppressed by the 
form factor, which is taken to be the one introduced by 
Helm [17]. Cross sections are normalized to a single nu- 
cleon by taking the assumed spin-independent enhance- 
ment oc A 2 on the different elements of the target into 
account. Limits on the cross section of a coherent WIMP- 
nucleon scattering process are calculated using the opti- 
mum acceptance region derived above and Yellin's Opti- 
mum Interval method [5]. Figure 3 compares the limit 



from the optimum acceptance region with the limit ob- 
tained from this dataset by the CRESST collaboration. 
There, the acceptance is defined to contain 90% of the 
tungsten recoils in the energy interval [10; 40] keV [4]. For 
comparison, some limits from other experiments are also 
shown. 
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FIG. 3: Exclusion plot on the WIMP-nucleon scattering cross 
section as function of WIMP mass, demonstrating the impact 
of the optimum acceptance region. The thick solid line (green) 
is the limit obtained with the optimum acceptance regions 
from the CRESST-II dataset discussed here as an example. 
The thick dashed line (light green) is the limit from the same 
dataset obtained by a one-sided 90% tungsten recoil accep- 
tance region in the energy interval [10;40]keV [4]. Limits 
from some other experiments are also shown: CRESST-I [18], 
CoGeNT [19], CDMS [7] and XENON10 [9]. 

At high WIMP masses, the optimum acceptance region 
gives a limit that is about 20% stronger than with the 
acceptance region previously used by the CRESST col- 
laboration. For low WIMP masses, the optimum accep- 
tance region results in drastically improved limits. For 
example, for a 10GeV/c 2 WIMP the improved limit is 
more than five orders of magnitude stronger than the 
one derived from a constant acceptance, that is to say, 
the optimum acceptance region allows to infer informa- 
tion about WIMPs in otherwise completely inaccessible 
regions of parameter space. 



III. COMBINING DATA FROM DIFFERENT 
DETECTORS 

The optimum acceptance region will be distinct for 
each detector, given its particular background and res- 
olution. Thus, the question arises for segmented-target 
experiments how to combine the individual subdetectors 
within the Maximum Gap or Optimum Interval meth- 
ods. Within these methods, the computed limit depends 
on integrals over the accepted signal spectrum dL a / dE 
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of the form 



rii := 



1+1 d£a 
dE 



dE 



(7) 



between two observed events in the acceptance region at 
energies Ei. To facilitate the combination of different 
detectors, the energy coordinate E is transformed to a 
new variable 77 according to the integral transform 



of the exposure and acceptance of each detector. There- 
fore, for each expected WIMP spectrum (and in par- 
ticular for each WIMP mass), individual detectors can 
now be joined to give the combined limit: The desired 
Maximum Gap or Optimum Interval method is simply 
applied to the energy interval [0, 1] considering the ob- 
served events from all detectors. The summed exposure 
is obtained by adding the individual M up. 



V(E) ■= 



dr a 

dE' 



dE' 



(8) 



IV. CONCLUSIONS 



which is a bijective transformation that leaves the cal- 
culated limit unchanged and occurs naturally within the 
Maximum Gap or Optimum Interval frameworks. The 
normalization constant M is chosen as the accepted- 
spectrum weighted exposure 



TV = 



d£a 

dE' 



dE'. 



(9) 



In this new energy variable 77, the accepted signal spec- 
trum is just a constant, which can be most easily seen 
manipulating differentials: 
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dE 



dr a 
dE 



dr a 

dE 



d 1 
dE AT 



1 



dTa 
dE 
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With this transformation, all possible differences of in- 
dividual detectors have been mapped into the interval 
77 = [0, 1] and the single number TV. Events of all detec- 
tors are distributed within 77 = [0, 1], and Af is a measure 



Previously, direkt dark matter searches would con- 
strain various WIMP models using one common accep- 
tance region. Here it was shown that by optimizing the 
acceptance region for each WIMP model, one can im- 
prove the sensitivity of an experiment by orders of mag- 
nitude. This has been demonstrated with recent data 
from CRESST-II as an example, where a drastically im- 
proved limit resulted in particular for low mass WIMPs. 
At the same time, the algorithm introduced here removes 
the ambiguity from defining the acceptance region in a 
rather ad hoc way. It was shown how to make full use of 
this optimization within the Maximum Gap or Optimum 
Interval frameworks to achieve a combined limit for in- 
dividual subdetectors of a segmcntcd-target experiment, 
or different experiments altogether. 

Acknowledgements 

I am thankful to my colleagues Jens Schmaler, Dieter 
Hauff and Franz Probst for useful discussions. 



[1] C. Amsler et al., The Review of Particle Physics, Physics 
Letters B 667 (2008), no. 1, available from http : //pdg . 
lbl.gov/, see in particular the Dark Matter review. 

[2] G. Jungman, M. Kamionkowski and K. Griest, Physics 
Reports 267 (1996), 195, arXiv :hep-ph/9506380. 

[3] R. J. Gaitskell, Annual Review of Nuclear and Particle 
Science 54 (2004), 315. 

[4] G. Angloher et al. (The CRESST Collaboration), As- 
troparticle Physics 31 (2009), 270, arXiv: 0809. 1829. 

[5] S. Yellin, Physical Review D 66 (2002), 032005, 
arXiv: physics/0203002, see also S. Yellin, 
arXiv: 0709. 2701 (2007). 

[6] G. Angloher et al. (The CRESST Collabora- 
tion), Astroparticle Physics 23 (2005), 325, 
arXiv : astro-ph/0408006. 

[7] Z. Ahmed et al. (The CDMS Collaboration), Physical 
Review Letters 102 (2009), 011301, arXiv: 0802. 3530. 

[8] V. Sanglard et al. (The EDELWEISS Collabo- 
ration), Physical Review D 71 (2005), 122002, 
arXiv : astro-ph/0503265. 

[9] J. Angle et al. (The XENON Collaboration), Physical 
Review Letters 100 (2008), 021303, arXiv: 0706. 0039. 
[10] G. J. Alner et al, Astroparticle Physics 28 (2007), 287, 



arXiv : astro-ph/0701858. 

[11] R. F. Lang et al., accepted for publication in Astroparti- 
cle Physics (2009), arXiv : 0905 .4282. 

[12] R. F. Lang et al., arXiv: 0910. 4414 (2009). 

[13] F. Donato, N. Fornengo and S. Scopel, Astroparticle 
Physics 9 (1998), 247, arXiv:hep-ph/9803295. 

[14] J. D. Lewin and P. F. Smith, Astroparticle Physics 
6 (1996), 87, see also Detector response correc- 
tions: correction and Spin factors - revised ta- 
bles, available from http : / /hepwww . rl . ac . uk/UKDMC/ 
pub/publications .html. 

[15] G. J. Feldman and R. D. Cousins, Physical Review D 57 
(1998), 3873, arXiv:physics/9711021. 

[16] R. F. Lang and W. Seidel, New Journal of Physics 11 
(2009), 105017, arXiv: 0906. 3290. 

[17] R. H. Helm, Physical Review 104 (1956), 1466. 

[18] G. Angloher et al. (The CRESST Collaboration), As- 
troparticle Physics 18 (2002), 43. 

[19] C. E. Aalseth et al. (The CoGeNT Collabora- 
tion), Physical Review Letters 101 (2008), 251301, 
arXiv : 0807 . 0879, see also erratum in Physical Review 
Letters 102 (2009), 109903(E). 



